WHAT IS CLAIMED IS: 



1 . A method for analyzing a computed tomography scan of a whole lung for lung 
nodules, said method comprising the steps of: 

(a) segmenting a first lung region and a second lung region from the 
computed tomography scan, the first lung region corresponding to lung parenchyma 
of the lung and the second lung region corresponding to at least one of a pleural 
surface of the lung and a surface defined by vessels within the lung; 

(b) generating an initial list of nodule candidates from the computed 
tomography scan within the first lung region, the list including at least a center 
location and an estimated size associated with each nodule candidate; 

(c) generating a subimage for each nodule candidate in the initial list; 

(d) selectively removing streaking artifacts from the subimage; and 

(e) filtering the nodule candidate identified on the initial list to eliminate 
false positives from the list. 

2. A method for analyzing a whole lung computed tomography scan as defined in 
Claim 1, wherein step (b) comprises the substeps of: 

(i) thresholding the first lung region; 

(ii) labeling high density voxels foreground voxels to identify nodule 
candidate regions; 

(iii) determining for each foreground voxel RMI; 

A 

(iv) selecting the local maximum R mi within a nodule candidate region; 

(v) determining a limited extent criterion for each foreground voxel which 

A 

corresponds to a /? Mil 
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(vi) generating the initial list of nodule candidates which satisfy the limited 

A 

extent criterion, the list including at least Nc and R mi associated with the 
corresponding foreground voxel. 

5 3. A method for analyzing a whole lung computed tomography scan as defined in 
Claim 1, wherein step (d) comprises the substeps of: 

(i) determining an amount of streaking artifact present in the sub-image; 

and 

10 (ii) filtering the streaking artifact out from the subimage when the amount 

of the streaking artifact present in the sub-image exceeds Tsar- 

4. A method for analyzing a whole lung computed tomography scan as defined in 
Claim 3, wherein the amount of streaking artifact present in the sub-image is 

1 5 calculated by a metric 

2 n m p 
Sm = S Z Z ( - KU + hk)f 
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5. A method for analyzing a whole lung computed tomography scan as defined in 
20 Claim 3, wherein the filtering is performed by a vertical median filter of size 1x3. 

6. A method for analyzing a whole lung computed tomography scan as defined in 
Claim 3, wherein Tsar is from about 20000 to about 80000. 

25 7. A method for analyzing a whole lung computed tomography scan as defined in 
Claim 1, wherein step (e) comprises the substeps of: 

(i) determining for each nodule candidate a fraction. Fa, a surface of the 
nodule candidate that is attached to other solid structures; and 
30 (ii) removing the nodule candidate from the list when the fraction exceeds 

Ta. 
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8. A method for analyzing a whole lung computed tomography scan as defined in 
Claim 1, wherein step (e) comprises the substeps of: 

5 (i) generating a cube wall about each nodule candidate; 

(ii) determining an intersection volume, Vnh corresponding to portions of 
the nodule region associated with the nodule candidate that intersect the cube 
wall; 

10 

(ii) removing the nodule candidate from the list when the firaction of the 
intersection volume, F„„ over the volume of the nodule candidate, V„, exceeds Tyy, 

9. A method for analyzing a whole lung computed tomography scan as defined in 
1 5 Claim 7, wherein step (e) comprises the substeps of: 

(i) generating a cube wall about each nodule candidate; 

(ii) determining an intersection volume, F„/, corresponding to portions of 
20 the nodule region associated with the nodule candidate that intersect the cube 

wall; 

(ii) removing the nodule candidate firom the list when the firaction of the 
intersection volume, Vnh over the volume of the nodule candidate, exceeds Tw. 

25 

10. A lung nodule detecting apparatus for analyzing a computed tomography scan 
of a whole lung for lung nodules, the lung nodule detecting apparatus comprising: 

a detecting unit configured to: 

30 

(a) segment a first lung region and a second lung region from the 
computed tomography scan, the first lung region corresponding to lung parenchyma 
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of the lung and the second lung region corresponding to at least one of a pleural 
surface of the lung and a surface defined by vessels within the lung; 

(b) generate an initial hst of nodule candidates from the computed 
5 tomography scan within the first lung region, the list including at least a center 

location and an estimated size associated with each nodule candidate; 

(c) generate a subimage for each nodule candidate in the initial list; 

10 (d) selectively remove streaking artifacts from the subimage; and 

(e) filter the nodule candidate identified on the initial list to eliminate false 
positives from the list. 

15 11. A lung nodule detecting apparatus as defined by Claim 10, wherein for step 
(b) said detecting unit is configured to: 

(i) threshold the first lung region; 

(ii) label high density voxels foreground voxels to identify nodule 
20 candidate regions; 

(iii) determine for each foreground voxel RMI; 

A 

(iv) select the local maximum R mi within a nodule candidate region; 

(v) determine a limited extent criterion for each foreground voxel which 

A 

corresponds toa R Mil 

25 (vi) generate the initial list of nodule candidates which satisfy the limited 

A 

extent criterion, the list including at least Nc and R mi associated with the 
corresponding foreground voxel. 

12, A lung nodule detecting apparatus as defined by Claim 10, wherein step (d) 
30 comprises the substeps of: 
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(i) determine an amount of streaking artifact present in the sub-image; and 

(ii) filter the streaking artifact out from the subimage when the amount of 
the streaking artifact present in the sub-image exceeds Tsar- 

13. A lung nodule detecting apparatus as defined by Claim 12, wherein the 
amount of streaking artifact present in the sub-image is calculated by a metric 



14. A lung nodule detecting apparatus as defined by Claim 12, further comprising 
a vertical median filter of size 1x3 to filter the streaking artifact. 

15. A lung nodule detecting apparatus as defined by Claim 12, wherein Tsar is 
from about 20000 to about 80000. 

16. A lung nodule detecting apparatus as defined by Claim 10, wherein step (e) 
comprises the substeps of: 

(i) determining for each nodule candidate a fraction. Fa, a surface of the 
nodule candidate that is attached to other solid structures; and 

(ii) removing the nodule candidate from the list when the fraction exceeds 



17. A lung nodule detecting apparatus as defined by Claim 10, wherein step (e) 
comprises the substeps of: 

(i) generating a cube wall about each nodule candidate; 

(ii) determining an intersection volume, Vnh corresponding to portions of 
the nodule region associated with the nodule candidate that intersect the cube 
wall; 




Ta. 
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(ii) removing the nodule candidate from the list when the fraction of the 
intersection volume, over the volume of the nodule candidate, F^, exceeds Tw. 

18. A lung nodule detecting apparatus as defined by Claim 16, wherein step (e) 
comprises the substeps of: 

(i) generating a cube wall about each nodule candidate; 

(ii) determining an intersection volume, V„i, corresponding to portions of 
the nodule region associated with the nodule candidate that intersect the cube 
wall; 

(ii) removing the nodule candidate from the list when the fraction of the 
intersection volume, F„i, over the volume of the nodule candidate, F^, exceeds Tyy. 

1 9. An article of manufacture for detecting lung nodules in a computed 
tomography scan of a whole lung, the article comprising: 

a machine readable medium containing one or more programs which when 
executed implement the steps of: 

(a) segmenting a first lung region and a second lung region from the 
computed tomography scan, the first lung region corresponding to lung parenchyma 
of the lung and the second lung region corresponding to at least one of a pleural 
surface of the lung and a surface defined by vessels within the lung; 

(b) generating an initial list of nodule candidates from the computed 
tomography scan within the first lung region, the list including at least a center 
location and an estimated size associated with each nodule candidate; 

(c) generating a subimage for each nodule candidate in the initial list; 
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selectively removing streaking artifacts from the subimage; and 



(e) filtering the nodule candidate identified on the initial list to eliminate 
false positives fi*om the list. 

5 

20. An article of manufacture for detecting lung nodules as defined by Claim 19, 
wherein step (b) comprises the substeps of: 

(i) thresholding the first lung region; 

(ii) labeling high density voxels foreground voxels to identify nodule 
candidate regions; 

(iii) determining for each foreground voxel RMI; 

A 

(iv) selecting the local maximum R mi within a nodule candidate region; 

(v) determining a limited extent criterion for each foreground voxel which 

A 

corresponds toa R Mil 

(vi) generating the initial list of nodule candidates which satisfy the limited 

A 

extent criterion, the list including at least Nc and R mi associated with the 
corresponding foreground voxel. 

20 21. An article of manufacture for detecting lung nodules as defined by Claim 19, 
wherein step (d) comprises the substeps of: 

(i) determining an amount of streaking artifact present in the sub-image; 

and 

25 (ii) filtering the streaking artifact out from the subimage when the amount 

of the streaking artifact present in the sub-image exceeds Tsar- 

22. An article of manufacture for detecting lung nodules as defined by Claim 21 , 
wherein the amount of streaking artifact present in the sub-image is calculated by a 
30 metric 
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23. An article of manufacture for detecting lung nodules as defined by Claim 21 , 
wherein the filtering is performed by a vertical median filter of size 1x3. 

24. An article of manufacture for detecting lung nodules as defined by Claim 21 , 
wherein T^ar is from about 20000 to about 80000. 



25. An article of manufacture for detecting lung nodules as defined by Claim 19, 
10 wherein step (e) comprises the substeps of; 

(i) determining for each nodule candidate a fraction, F^,, a surface of the 
nodule candidate that is attached to other solid structures; and 

(ii) removing the nodule candidate from the list when the fraction exceeds 

15 r«. 

26. An article of manufacture for detecting lung nodules as defined by Claim 19, 
wherein step (e) comprises the substeps of: 

20 (i) generating a cube wall about each nodule candidate; 

(ii) determining an intersection volume, Vnh corresponding to portions of 
the nodule region associated with the nodule candidate that intersect the cube 
wall; 

25 

(ii) removing the nodule candidate from the list when the fraction of the 
intersection volume, Vnu over the volume of the nodule candidate, Vn, exceeds 
T ^. 

30 27. An article of manufacture for detecting lung nodules as defined by Claim 25, 
wherein step (e) comprises the substeps of: 
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(i) generating a cube wall about each nodule candidate; 

(ii) determining an intersection volume, F„„ corresponding to portions of 
the nodule region associated with the nodule candidate that intersect the cube 

5 wall; 

(ii) removing the nodule candidate from the list when the fraction of the 
intersection volume, V„ij over the volume of the nodule candidate, V„, exceeds T^. 

28. A method for correlating a segmentation of 3-d images of a pulmonary nodule 
from a high-resolution computed tomography (CT) scans, the images being in a 
floating point pixel-format associated with a 6-dimensional parameter space and 

10 including a first image (im\) obtained at time-1 and a second image {imz) obtained at 
time-2, the method comprising the steps of: 

(a) selecting a first region-of-interest (ROIi) for the nodule in the first 
image (im\); 

(b) selecting a second region-of-interest (ROI2) for the nodule in the 
1 5 second image (/mi); 

(c) registering the second region-of-interest (ROI2) to the first region-of- 
interest (ROIi) to obtain a transformed second region-of-interest (ROIir); 

(d) separately segmenting both the nodule in the first region-of-interest 
(ROIi) and the transformed second region-of-interest (ROI2/); and 

20 (e) adjusting the first segmented nodule (5i) and the second segmented 

nodule (Si)- 

29. A method for correlating a segmentation of 3-d images as defined in Claim 28, 
wherein the first region-of-interest (ROIi) is cubic and is selected to be about three 

25 times the size of the diameter of the nodule. 
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30. A method for correlating a segmentation of 3-d images as defined in Claim 28, 
wherein the second region-of-interest (ROI2) is cubic and is selected to be about three 
times the size of the diameter of the nodule. 

5 31. A method for correlating a segmentation of 3-d images as defined in Claim 28, 
wherein step (d) includes at least one of the following substeps of: 

(i) gray-level thresholding; 

(ii) morphological filtering for vessel removal: and 

(iii) plane clipping for separarting a pleural wall. 

10 

32. A method for correlating a segmentation of 3-d images as defined in Claim 3 1 , 
wherein the gray-level thresholding is performed at an adaptive threshold level. 

33. A method for correlating a segmentation of 3-d images as defined in Claim 32, 
15 wherein the adaptive threshold level is selected for each region-of-interest (ROIi and 

ROI2) by: 

determining a peak parenchyma value, Vpi 
determining a peak nodule value, v^,; 

calculating the adaptive threshold level as a midpoint between the peak 
20 parenchyma value, v^, and the peak nodule value, v„. 

34. A method for correlating a segmentation of 3-d images as defined in Claim 33, 
further comprising the step of calculating an intensity histogram, H(x) for 
determining the peak parenchyma value, Vp, and the peak nodule value, v„. 

25 

35. A method for correlating a segmentation of 3-d images as defined in Claim 34, 
wherein the intensity histogram, H(x)^ is calculated between about -1024 HU and 
about 476 HU with a bin size of about 1 . 
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36. A method for correlating a segmentation of 3-d images as defined in Claim 34, 
further comprising the substep of: 

filtering the intensity histogram, H{x), with a gaussian with standard 
5 deviation of about 25 HU prior to determining the peak parenchyma value, Vp, and 
the peak nodule value, v„. 

37. A method for correlating a segmentation of 3-d images as defined in Claim 34, 
wherein the intensity histogram, H{x), is searched between about -1024 HU and 

10 about -680 for the peak parenchyma value, Vp. 

38. A method for correlating a segmentation of 3-d images as defined in Claim 34, 
wherein the intensity histogram, //(jc), is searched between about -227 HU and 
about -173 for the peak nodule value, v„. 

15 

39. A method for correlating a segmentation of 3-d images as defined in Claim 28, 
wherein registering the second region-of-interest (ROI2) to the first region-of-interest 
(ROIi) to obtain a transformed second region-of-interest (ROI2O comprises the 
substeps of: 

20 (a) calculating initial rigid-body transformation parameters for a rigid- 

body transformation on the second region-of-interest (ROI2); 

(b) determining tfie optimum rigid-body transformation parameters by 
calculating a registration metric between the first region-of-interest (ROIi) and the 
rigid-body transformation on the second region-of-interest (ROI2); and 

25 (c) generating a registered image firom the optimum rigid-body 

transformation parameters. 
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40. A method for correlating a segmentation of 3-d images as defined in Claim 39, 
wherein the registration metric is calculated by 

transforming the second region-of-interest (ROI2) with the initial rigid- 
body transformation parameters to obtain a transformed second region-of-interest 
(ROI2,); 

calculating the registration metric as a mean-squared-difference (MSD) 
between the transformed second regibn-of-interest (ROI2/) and the first region-of- 
interest (ROIi); and 

searching for the minimum mean-squared-difiference (MSD) in the 6- 
dimensional parameter space. 
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4 L A method for correlating a segmentation of 3-d images as defined in Claim 40, 
wherein the transforming of the second region-of-interest (ROI2) to obtain the 
transformed second region-of-interest (ROl2r) is a mapping of a point v in 3-d space 
to a point v' in transformed space defined by : 



L^zJ 



wherein R^, Ry, and R^ are rotation matrices defined as: 
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R. = 



R.. = 



10 0 
0 cos(r,) -sin(r,) 
0 sin(rj cos(rJ 



cos(r^) 0 sin(r^) 

0 1 0 
-sin(r^) 0 cos(r^) 
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R. = 



"cos(rJ 
sin(rj 
0 



-sm(rj O' 
cos(r^ ) 0 
0 1 



42. A method for correlating a segmentation of 3-d images as defined in Claim 39, 
wherein said initial rigid-body transformation parameters include six parameters 
(tx,ty,tz,rx,ry,rz) respectively defined as translation in x, translation in y, translation in 
z, rotation about the x-axis, rotation about the y-axis, and rotation about the z-axis; 

wherein the initial rotation parameters (rx,ry,rz) are all set to zero; and the 
initial translation parameters (tx,ty,tz,) are set so that the nodule in the first region-of- 
interest (ROIi) overlaps the nodule in the second region-of-interest (ROI2)) during the 
initial calculation of the registration metric. 



43. A method for correlating a segmentation of 3-d images as defined in Claim 40, 
wherein the mean-squared-difference (MSD) is gaussian weighted. 



44. A method for correlating a segmentation of 3-d images as defined in Claim 31, 
wherein: 

a first thresholded image (ri)and a second thresholded image (T2) are defined 
by gray-level thresholding prior to vessel removal and separating the pleural wall; 
and 

step (e) is performed by comparing the segmented nodules and the thresholded 

images. 



45. A method for correlating a segmentation of 3-d images as defined in Claim 44, 
wherein an active pixel is marked as one of: 

a repeat nodule pixel; 

a nodule growth pixel; 
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a nodule atrophy pixel; and 

a nodule missegmentation pixel. 

46. A method for correlating a segmentation of 3-d images as defined in Claim 44, 
wherein a foreground pixel in the first segmented nodule (^i) is marked as a repeated 
nodule pixel from the first region-of-interest (ROIi) to the transformed second region- 
of-interest (ROI2/) when the corresponding pixel in second segmented nodule (^2) and 
the corresponding pixel in second thresholded image (T2) are both foreground 

47. A method for correlating a segmentation of 3-d images as defined in Claim 45, 
wherein a foreground pixel in the first segmented nodule (Si) is marked as a nodule 
atrophy pixel when the corresponding pixel in second segmented nodule (82) is 
background and the corresponding pixel in second thresholded image (T2) is 
background. 

48. A method for correlating a segmentation of 3-d images as defined in Claim 45, 
wherein a foreground pixel in the first segmented nodule (S\) is marked as a 
missegmented pixel in the first region-of-interest (ROIi) when the corresponding 
pixel in second segmented nodule (S2) is background and the corresponding pixel in 
second thresholded image (72) is foreground. 

49. A method for correlating a segmentation of 3-d images as defined in Claim 45, 
wherein a foreground pixel in the second segmented nodule (52) is marked as a 
repeated nodule pixel from the first region-of-interest (ROIi) to the transformed 
second region-of-interest (ROI2/) when the corresponding pixel in first segmented 
nodule (JS\) and the corresponding pixel in first thresholded image (T\) are both 
foreground. 
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50. A method for correlating a segmentation of 3-d images as defined in Claim 45, 
wherein a foreground pixel in the second segmented nodule (S2) is marked as a nodule 
growth pixel when the corresponding pixel in first segmented nodule (Si) is 
background and the corresponding pixel in first thresholded image (T\) is background. 

5 

51 . A method for correlating a segmentation of 3-d images as defined in Claim 45, 
wherein a foreground pixel in the second segmented nodule (S2) is marked as a 
missegmented pixel in the transformed second region-of-interest (ROI2/) when the 
corresponding pixel in first segmented nodule (^i) is background and the 

10 corresponding pixel in first thresholded image (T\) is foreground. 

52. A registration apparatus for correlating a segmentation of 3-d images of a 
pulmonary nodule firom a high-resolution computed tomography (CT) scans, the 
images being in a floating point pixel-format associated with a 6-dimensional 

15 parameter space and including a first image (im\) obtained at time-1 and a second 
image (zmi) obtained at time-2, the registration apparatus comprising: 

a registration unit configured to: 

(a) select a first region-of-interest (ROIi) for the nodule in the first image 

(irnx); 

20 (b) select a second region-of-interest (ROI2) for the nodule in the second 

image (1/W2); 

(c) register the second region-of-interest (ROI2) to the first region-of- 
interest (ROIi) to obtain a transformed second region-of-interest (ROI2/); 

(d) separately segment both the nodule in the first region-of-interest 
25 (ROIi) and the transformed second region-of-interest (ROIir); and 

(e) adjust the first segmented nodule and the second segmented 
nodule (82)- 
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53. A registration apparatus for correlating a segmentation of 3-d images as 
defined in Claim 52, wherein the first region-of-interest (ROIi) is cubic and is 
selected to be about three times the size of the diameter of the nodule. 



5 54. A registration apparatus for correlating a segmentation of 3-d images as 
defined in Claim 52, wherein the second region-of-interest (ROI2) is cubic and is 
selected to be about three times the size of the diameter of the nodule. 

55. A registration apparatus for correlating a segmentation of 3-d images as 
10 defined in Claim 52, wherein for step (d) said registration unit is configured to: 

(i) gray-level threshold; 

(ii) morphological filter to remove vessels; and 

(iii) plane clip to separate a pleural wall. 

15 56. A registration apparatus for correlating a segmentation of 3-d images as 
defined in Claim 55, wherein the gray-level threshold is performed at an adaptive 
threshold level. 

57. A registration apparatus for correlating a segmentation of 3-d images as 

20 defined in Claim 56, wherein the adaptive threshold level is selected for each region- 
of-interest (ROIi and ROI2) by: 

determining a peak parenchyma value, Vp\ 

determining a peak nodule value, v„; 

calculating the adaptive threshold level as a midpoint between the peak 
25 parenchyma value, Vp, and the peak nodule value, Vn- 

58. A registration apparatus for correlating a segmentation of 3-d images as 

defined in Claim 57, wherein said registration unit is configured to calculate an 
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intensity histogram, H(x) for determining the peak parenchyma value, Vp, and the 
peak nodule value, v„. 

59. A registration apparatus for correlating a segmentation of 3-d images as 

5 defined in Claim 58, wherein the intensity histogram, H(x), is calculated between 
about -1024 HU and about 476 HU with a bin size of about 1 . 

60. A registration apparatus for correlating a segmentation of 3-d images as 
defined in Claim 58, wherein said registration unit is configured to filter the intensity 

10 histogram, H(x), with a gaussian with standard deviation of about 25 HU prior to 
determining the peak parenchyma value, v^, and the peak nodule value, v„. 

61 . A registration apparatus for correlating a segmentation of 3-d images as 
defined in Claim 58, wherein the intensity histogram, H(x), is searched between 

1 5 about -1 024 HU and about -680 for the peak parenchyma value, Vp. 

62. A registration apparatus for correlating a segmentation of 3-d images as 
defined in Claim 58, wherein the intensity histogram, H(x), is searched between 
about -227 HU and about -173 for the peak nodule value, v„. 

20 

63. A registration apparatus for correlating a segmentation of 3-d images as 
defined in Claim 52, wherein said registration unit is configured to: 

(a) calculate initial rigid-body transformation parameters for a rigid-body 
transformation on the second region-of-interest (ROI2); 

25 (b) determine the optimum rigid-body transformation parameters by 

calculating a registration metric between the first region-of-interest (ROIi) and the 
rigid-body transformation on the second region-of-interest (ROI2); and 
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(c) generate a registered image from the optimum rigid-body 
transformation parameters. 



64. A registration apparatus for correlating a segmentation of 3-d images as 
defined in Claim 63, wherein the registration metric is calculated by 

transforming flie second region-of-interest (ROI2) with the initial rigid- 
body transformation parameters to obtain a transformed second region-of-interest 
(ROI2.)); 

calculating the registration metric as a mean-squared-difference (MSD) 
between the transformed second region-of-interest (ROI2O and the first region-of- 
interest (ROIi); and 

searching for the minimum mean-squared-difference (MSD) in the 6- 
dimensional parameter space. 



65. A registration apparatus for correlating a segmentation of 3-d images as 
defined in Claim 64, wherein the transforming of the second region-of-interest (ROI2) 
to obtain the transformed second region-of-interest (ROI2O is a mapping of a point v 
in 3-d space to a point v* in transformed space defined by : 

wherein R^^^ Ry, and R^ are rotation matrices defined as: 



V=R^RyR^v^ 



R. = 



10 0 
0 cos(rJ -sin(r,) 
0 sin(r^) cos(r^) 
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cos(r^) 
0 

-sin(r^) 



0 siriCr^,) 

1 0 

0 cos(r^) 



cos(rJ 
sin(rj 
0 



-sin(rj 0 
cos(r^ ) 0 
0 1 



66. A registration apparatus for correlating a segmentation of 3-d images as 
defined in Claim 63, wherein said initial rigid-body transformation parameters include 
six parameters (tx,tyjz,rx,ry,rz) respectively defined as translation in x, translation in 
y, translation in z, rotation about the x-axis, rotation about the y-axis, and rotation 
about the z-axis; 

wherein the initial rotation parameters (rx^ry^rz) are all set to zero; and the 
initial translation parameters (/jc,fy,fe,) are set so that the nodule in the first region-of- 
interest (ROIi) overlaps the nodule in the second region-of-interest (ROI2)) during the 
initial calculation of the registration metric. 

67. A registration apparatus for correlating a segmentation of 3-d images as 
defined in Claim 64, wherein the mean-squared-difference (MSD) is gaussian 
weighted. 

68. A registration apparatus for correlating a segmentation of 3-d images as 
defined in Claim 55, wherein: 

a first thresholded image (ri)and a second thresholded image ( J2) are defined 
by gray-level thresholding prior to vessel removal and separating the pleural wall; 
and 

step (e) is performed by comparing the segmented nodules and the thresholded 
images. 
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69. A registration apparatus for correlating a segmentation of 3-d images as 
defined in Claim 68, wherein an active pixel is marked as one of: 

a repeat nodule pixel; 

a nodule growth pixel; 

a nodule atrophy pixel; and 

a nodule missegmentation pixel. 

70. A registration apparatus for correlating a segmentation of 3-d images as 
defined in Claim 69, wherein a foreground pixel in the first segmented nodule is 
marked as a repeated nodule pixel firom the first region-of-interest (ROIi) to the 
transformed second region-of-interest (ROI2/) when the corresponding pixel in second 
segmented nodule (S2) and the corresponding pixel in second thresholded image {T-i) 
are both foreground 

71 . A registration apparatus for correlating a segmentation of 3-d images as 
defined in Claim 69, wherein a foreground pixel in the first segmented nodule is 
marked as a nodule atrophy pixel when the corresponding pixel in second segmented 
nodule {Si) is background and the corresponding pixel in second thresholded image 
(72) is background. 

72. A registration apparatus for correlating a segmentation of 3-d images as 
defined in Claim 69, wherein a foreground pixel in the first segmented nodule is 
marked as a missegmented pixel in the first region-of-interest (ROIi) when the 
corresponding pixel in second segmented nodule (^2) is background and the 
corresponding pixel in second thresholded image {Ti) is foreground. 
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73. A registration apparatus for correlating a segmentation of 3-d images as 
defined in Claim 69, wherein a foreground pixel in the second segmented nodule (S2) 
is marked as a repeated nodule pixel from the first region-of-interest (ROIi) to the 
transformed second region-of-interest (ROI2/) when the corresponding pixel in first 

5 segmented nodule (5i) and the corresponding pixel in first thresholded image (T\) are 
both foreground. 

74. A registration apparatus for correlating a segmentation of 3-d images as 
defined in Claim 69, wherein a foreground pixel in the second segmented nodule (S2) 

10 is marked as a nodule growth pixel when the corresponding pixel in first segmented 
nodule (S\) is background and the corresponding pixel in first thresholded image (Ji) 
is background. 

75. A registration apparatus for correlating a segmentation of 3-d images as 

1 5 defined in Claim 69, wherein a foreground pixel in the second segmented nodule (S2) 
is marked as a missegmented pixel in the transformed second region-of-interest 
(ROI2O when the corresponding pixel in first segmented nodule (S\) is background 
and the corresponding pixel in first thresholded image (T\) is foreground. 

20 76. An article of manufacture for correlating a segmentation of 3-d images of a 
pulmonary nodule from a high-resolution computed tomography (CT) scans, the 
images being in a floating point pixel-format associated with a 6-dimensional 
parameter space and including a first image (i/wi) obtained at time-1 and a second 
image (//W2) obtained at time-2, the article comprising: 

25 a machine readable medium containing one or more programs which when 

executed implement the steps of: 

(a) selecting a first region-of-interest (ROIi) for^the nodule in the first 
image (zmi); 

(b) selecting a second region-of-interest (ROI2) for the nodule in the 
30 second image (//W2); 
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(c) registering the second region-of-interest (ROI2) to the first region-of- 
interest (ROIi) to obtain a transformed second region-of-interest (ROIar); 

(d) separately segmenting both the nodule in the first region-of-interest 
(ROIi) and the transformed second region-of-interest (ROl2r); and 

(e) adjusting the first segmented nodule (S\) and the second segmented 
nodule {S2). 

77. An article of manufacture for correlating a segmentation of 3-d images as 
defined in Claim 76, wherein the first region-of-interest (ROIi) is cubic and is 
selected to be about three times the size of the diameter of the nodule. 

78. An article of manufacture for correlating a segmentation of 3-d images as 
defined in Claim 76, wherein the second region-of-interest (ROI2) is cubic and is 
selected to be about three times the size of the diameter of the nodule. 

79. An article of manufacture for correlating a segmentation of 3-d images as 
defined in Claim 76, wherein step (d) includes at least one of the following substeps 
of: 

(i) gray-level thresholding; 

(ii) morphological filtering for vessel removal: and 

(iii) plane clipping for separarting a pleural wall. 

80. An article of manufacture for correlating a segmentation of 3-d images as 
defined in Claim 79, wherein the gray-level thresholding is performed at an adaptive 
threshold level. 

81 . An article of manufacture for correlating a segmentation of 3-d images as 
defined in Claim 80, wherein the adaptive threshold level is selected for each region- 
of-interest (ROIi and ROI2) by: 
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determining a peak parenchyma value, Vp\ 
determining a peak nodule value, v„; 

calculating the adaptive threshold level as a midpoint between the peak 
parenchyma value, v^, and the peak nodule value, v„. 

5 

82. An article of manufacture for correlating a segmentation of 3-d images as 
defined in Claim 81, fiirther comprising the step of calculating an intensity histogram, 
H{x) for determining the peak parenchyma value, v^, and the peak nodule value, v„. 

10 83. An article of manufacture for correlating a segmentation of 3-d images as 
defined in Claim 82, wherein the intensity histogram, H{pc)^ is calculated between 
about -1024 HU and about 476 HU with a bin size of about 1 . 

84. An article of manufacture for correlating a segmentation of 3-d images as 
15 defined in Claim 82, further comprising the substep of: 

filtering the intensity histogram, H{x)^ with a gaussian with standard 
deviation of about 25 HU prior to determining the peak parenchyma value, Vp, and 
the peak nodule value, v„. 

20 85. An article of manufacture for correlating a segmentation of 3-d images as 
defined in Claim 82, wherein the intensity histogram, //(x), is searched between 
about -1024 HU and about -680 for the peak parenchyma value, Vp. 

86. An article of manufacture for correlating a segmentation of 3-d images as 
25 defined in Claim 82, wherein the intensity histogram, H{x\ is searched between 
about -227 HU and about -173 for the peak nodule value, v„. 
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87. An article of manufacture for correlating a segmentation of 3-d images as 
defined in Claim 76, wherein registering the second region-of-interest (ROI2) to the 
first region-of-interest (ROIi) to obtain a transformed second region-of-interest 
(ROl2r) comprises the substeps of: 

(a) calculating initial rigid-body transformation parameters for a rigid- 
body transformation on the second region-of-interest (ROI2); 

(b) determining the optimum rigid-body transformation parameters by 
calculating a registration metric between the first region-of-interest (ROIi) and the 
rigid-body transformation on the second region-of-interest (ROI2); and 

(c) generating a registered image from the optimum rigid-body 
transformation parameters. 

88. An article of manufacture for correlating a segmentation of 3-d images as 
defined in Claim 87, wherein the registration metric is calculated by 

transforming the second region-of-interest (ROI2) with the initial rigid- 
body transformation parameters to obtain a transformed second region-of-interest 

(ROI2O); 

calculating the registration metric as a mean-squared-difference (MSD) 
between the transformed second region-of-interest (ROI2/) and the first region-of- 
interest (ROIi); and 

searching for the minimum mean-squared-difference (MSD) in the 6- 
dimensional parameter space. 

89. An article of manufacture for correlating a segmentation of 3-d images as 
defined in Claim 88, wherein the transforming of the second region-of-interest (ROI2) 
to obtain the transformed second region-of-interest (ROI2/) is a mapping of a point v 
in 3-d space to a point v' in transformed space defined by : 
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wherein Ry^ and R^ are rotation matrices defined as: 



R. = 



10 0 
0 cos(r,) -sin(r^) 
0 sin(r^) cos(r^) 



R. = 



cos(r^) 0 sin(r^) 

0 1 0 
-sin(r^) 0 cos(r^)j 



R, = 



cos(rJ 
sin(rj 
0 



-sin(rj 0 
cos(rj ) 0 
0 1 



90. An article of manufacture for correlating a segmentation of 3-d images as 
defined in Claim 87, wherein said initial rigid-body transformation parameters include 
six parameters (tXyty,tz,rx,ry,rz) respectively defined as translation in x, translation in 
y, translation in z, rotation about the x-axis, rotation about the y-axis, and rotation 
about the z-axis; 

wherein the initial rotation parameters {rx^ry^rz) are all set to zero; and the 
initial translation parameters (tx,ty,tZy) are set so that the nodule in the first region-of- 
interest (ROIi) overlaps the nodule in the second region-of-interest (ROI2)) during the 
initial calculation of the registration metric. 



91 . An article of manufacture for correlating a segmentation of 3-d images as 
defined in Claim 88, wherein the mean-squared-difference (MSD) is gaussian 
weighted. 



64 



92. An article of manufacture for correlating a segmentation of 3-d images as 
defined in Claim 79, wherein: 

a first thresholded image (ri)and a second thresholded image (T2) are defined 
by gray-level thresholding prior to vessel removal and separating the pleural wall; 
and 

step (e) is performed by comparing the segmented nodules and the thresholded 

images. 

93. An article of manufacture for correlating a segmentation of 3-d images as 
defined in Claim 92, wherein an active pixel is marked as one of: 

a repeat nodule pixel; 

a nodule growth pixel; 

a nodule atrophy pixel; and 

a nodule missegmentation pixel. 

94. An article of manufacture for correlating a segmentation of 3-d images as 
defined in Claim 93, wherein a foreground pixel in the first segmented nodule (Si) is 
marked as a repeated nodule pixel from the first region-of-interest (ROIi) to the 
transformed second region-of-interest (ROI2O when the corresponding pixel in second 
segmented nodule (^2) and the corresponding pixel in second thresholded image (72) 
are both foreground 

95. An article of manufacture for correlating a segmentation of 3-d images as 
defined in Claim 93, wherein a foreground pixel in the first segmented nodule is 
marked as a nodule atrophy pixel when the corresponding pixel in second segmented 
nodule (^2) is background and the corresponding pixel in second thresholded image 
(72) is background. 
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96. An article of manufacture for correlating a segmentation of 3-d images as 
defined in Claim 93, wherein a foreground pixel in the first segmented nodule (JSi) is 
marked as a missegmented pixel in the first region-of-interest (ROIi) when the 

5 corresponding pixel in second segmented nodule (^2) is background and the 
corresponding pixel in second thresholded image (T2) is foreground. 

97. An article of manufacture for correlating a segmentation of 3-d images as 
defined in Claim 93, wherein a foreground pixel in the second segmented nodule (S2) 

10 is marked as a repeated nodule pixel fi*om the first region-of-interest (ROIi) to the 
transformed second region-of-interest (ROI2/) when the corresponding pixel in first 
segmented nodule (Si) and the corresponding pixel in first thresholded image (T\) are 
both foreground. 

15 98. An article of manufacture for correlating a segmentation of 3-d images as 

defined in Claim 93, wherein a foreground pixel in the second segmented nodule (S2) 
is marked as a nodule growth pixel when the corresponding pixel in first segmented 
nodule (S\) is background and the corresponding pixel in first thresholded image (T\) 
is background. 

20 

99. An article of manufacture for correlating a segmentation of 3-d images as 
defined in Claim 93, wherein a foreground pixel in the second segmented nodule (S2) 
is marked as a missegmented pixel in the transformed second region-of-interest 
(ROl2r) when the corresponding pixel in first segmented nodule (S\) is background 
25 and the corresponding pixel in first thresholded image (T\) is foreground. 
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